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Dielectric Relaxation for Spherical Molecules in a 
Crystalline Field: Theory for Two Simple Models 

John D. Hoffman and Benjamin M. Axilrod 

A theory of dielectric relaxation is given for two simple models where the dipoles may 
occupy orientational sites arranged in three dimensions. In the first model, the dipole may 
point toward the corners of a regular tetrahedron, and in the second, toward the faces of a 
cube. Both models are simplified to the extent that only one of the sites can differ in energy 
from all the others. 

The calculations show that these models lead to a discrete set of dielectric relaxation 
times when the sites are not equivalent. The effect of the set of relaxation times on the 
characteristics of the dielectric relaxation spectrum is examined. The relaxation times are 
always sufficiently close together so that the dielectric loss is confined to a single region. 
However, the separation of the relaxation times is sufficient to cause the dielectric loss peak 
to be considerably broader than that predicted by the Debye equations. In addition, it is 
found that the loss maximum will frequently be asymmetrical. The variation of the di- 
electric constant with temperature, and the breadth and asymmetry of the loss peaks are 
discussed. 



1. Introduction 

Our basic viewpoint is that an anisotropy in the 
crystalline field hindering the rotation of a dipole in 
a lattice can be the source of a set of discrete dielectric 
relaxation times. This viewpoint has already been 
explored in some detail for a single-axis rotator. 1 2 
The main objective of this paper is to show how one 
type of calculation used with the single-axis rotator 
can be extended to include two simple models where 
each dipole in the lattice, by turning about a lattice 
point, can occupy orientational sites arranged in 
three dimensions. A model permitting such ori- 
entations in space may be a reasonable representation 
for certain molecular crystals composed of polar 
molecules that are approximately spherical. 

The theory presented here and in the papers men- 
tioned above may be regarded as an extension of 



i J. D. Hoffman and H. G. Pfeiffer, J. Chem. Phys. 22, 132 (1954). 
2 J. D. Hoffman, ibid, (in press). 



Debye's two-position model. 3 The two-position 
model gives only a single l 3 4 relaxation time ; how- 
ever, a set of relaxation times frequently appears for 
a single-axis rotator when more than two sites are 
involved. 1 2 We have also modified Debye's theory 
to the extent that transition probabilities of the 
elementary process for dipole reorientation have been 
defined along the lines proposed by Kauzmann. 5 

2. Tetrahedral Model With One Stable Site 

2.1. Model 

In this model the dipole possesses four equilibrium 
orientations (sites) designated 1 through 4 which 
are arranged in a tetrahedral manner (fig. 1). The 
permanent dipole moment has the magnitude fx. 

3 P. Debye, Polar Molecules, p. 104 (Dover Publications, New York, 1945). 

4 H. Frohlich, Theory of Dielectrics, p. 83 (Oxford University Press, 1949). 
s W. Kauzmann, Rev. Mod. Phys. 14, 12 (1942). 





0° 109.5° 
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Figure 1. Elementary processes for the tetrahedral model. 

Left diagram: Heavy straight arrows indicate equilibrium dipole orientations, light curved arrows the elementary 

processes k ( >) and k r ( — ►)• Site 1 is the stable orientation. Right diagram: Details of the barrier system for site 1 

relative to sites 2, 3, and 4. The barrier between sites 2, 3, and 4 (not shown) is W, and the corresponding transition 
probability is k'. 
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Sites 2, 3, and 4 are equivalent; site 1 is unique, and 
may be more stable than the others. 

The dipole can turn in a single jump (elementary 
process) from site 1 to either site 2, 3, or 4 by sur- 
mounting a local free energy of activation barrier 
W-f-V. W and V are shown in figure 1 (right 
diagram). The barrier system is considered to be 
the result of the crystalline field imposed on the 
dipole by its neighbors in the lattice. It is assumed 
in this simplified treatment that there is no correla- 
tion with the motion of neighboring dipoles. 

The probability that a dipole will jump from site 
1 to either site 2, 3, or 4 is given by 



k==B e -(W+V)lkT 



(1) 



where B is a frequency factor, k is Boltzmann's con- 
stant, and T the absolute temperature. The dipole 
must surmount an activation barrier W in order to 
turn from either site 2, 3, or 4 back to site 1, and the 
corresponding probability is 



!c'=B e- w l kT . 



(2) 



The barrier between any pair of the sites 2, 3, and 4 
is also W so that the probability of the reorientations 
2^3, 3^4, and 2^4 is k'. ft is readily seen that 
when F>0, so that site 1 is relatively stable, we 
have k f ^>k. The assigned transition probabilities 
for reorientations between the various sites are 
shown in figure 1 (left diagram.) 

2.2. Rate Equations 

The molecular relaxation times associated with the 
model are calculated by constructing and solving 
the differential equations that describe the rate at 
which dipoles enter and leave each site in terms of 
the elementary process transition probabilities given 
by eq (1) and (2). The number of dipoles leaving a 
site i and entering another designated j is the product 
of the probability of the transition i-*j and the pop- 
ulation, Nf, of site i. In our problem each k repre- 
sents the probability that a dipole will leave a site 
in a specified direction and turn to an adjacent site 
(single-jump hypothesis, S-J). 1 ' 2 

Consider the net rate at which dipoles enter and 
leave site 1. According to our hypothesis the num- 
ber per second leaving site 1 and entering site 2 is 
kNi) the total number leaving to enter sites 2, 3, 
and 4 is 3kN x . Similarly, the number returning 
from site 2 to site 1 is k'N 2 . The corresponding 
numbers leaving sites 3 and 4 to enter site 1 are 
k'Nz and k'N±. Hence the net rate for site 1 is 
given by dA\/dt= — (number leaving) + (number en- 
termg) = —3kN l +k / (N 2 +N s +N 4 ) . Similar reason- 
ing leads to the differential equations for the net 
rates for the other sites. Thus we find 

dN 1 /dt=-3kN 1 +k'(N 2 +N s +N 4 ) * 
dN 2 /dt=kN 1 -3k r N 2 +k'(N s +N.) 
dN 3 /dt=kN 1 +k , N 2 -3k / N s +k f N, 
dNJdt=kN 1 +k'(N 2 +N s )-3k'N i .^ 



The solutions of eq (3) will provide the molecular 
(intrinsic) relaxation times, Tp, where fi is the mode 
of decay, and will also give considerable information 
regarding the nature of the modes of decay. 

2.3. Relaxation Times and Modes of Decay 

Solutions of the equations of the type illustrated 
in eq (3) are of the form 1,2 - 6 







fl, 



(4) 



where i refers to a particular site, 12 is the number of 
sites, p is the mode of decay, and \[/p is a decay func- 
tion of the form 



fe= 



pi ' (A, *') t 



The set of discrete relaxation times is given by 



(5) 



(6) 



The solutions of eq (3) are obtained by forming the 
characteristic determinant of the differential equa- 
tions, and solving for eigenvalues of the operator 
D=d/dt. These eigenvalues, which are always found 
to be negative or zero, are the fp(k,k'). One of the 
fp(k,k f ) is always zero, and this mode is designated 
by 0=1; the corresponding C a are the equilibrium 
values of the N t . The relationships between the 
Ctfi for each mode where /3^2 are determined by sub- 
stituting each solution Ci$$ into eq (3). On carry- 
ing out the above steps we find that the solutions of 
eq (3) may be written 7 



Mode 2 Mode 3 

N k = C n + 3<? 22 tA 2 + 



n 2 =c 21 -- c*h + Ctfh 



M=C 4I -- r^ 



C 4 Z^3 



(7) 



where \l/ 2 = ex-p[—(k'-\-3k)f\ and ^ 3 =exp(— 4:k't). 
Also, we may note that C 2 i = C 3 i = C4i = C n exp(— V/ 
kT), and that C 23 + C p +C 4Z =0. 

Equation (7) describes the rate at which equilib- 
rium is established after the abrupt removal of a dis- 
turbance at t=0. In the present case the displace- 
ment from equilibrium at tS will be assumed to be 
the result of a static electric field. The quantities 
Ca and C i3 in eq (7) represent the (small) number of 
dipoles that are displaced from their equilibrium 
orientations by the field. 

If we form the vectors C i2 where each C i2 has the 
direction of a dipole in the it\\ site and the magnitude 



6 Ince, Ordinary Differential Equations, p. 144 (Dover Publications, New York, 
1944). 

7 The eigenvalue of the characteristic determinant D = — 4fc' occurs twice. In 
the case of such doubly degenerate roots we actually have Ni— Ci\-\-Civi/2-\-dypz 
-\-tCi\\pz. However, it is found by substituting this solution into the original 
differential equations that the Ca are identically zero. 
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Ci2, the sum of the C i2 is a measure of the polarization 
associated with mode 2. A similar relation applies 
to mode 3. Thus, eq (7) indicates that, owing to the 
return of the displaced dipoles to their equilibrium 
orientations, the polarization will decay with time 
according to 

P(t)=P 2 e-«' +mt +P 3 e- 4k ' t . (8) 

P 2 and P 3 are the polarizations associated with modes 
2 and 3 at /=0. The actual values of P> 2 and P 3 are 
calculated in a following section. Since the dielectric 
relaxation time is defined as the time the polarization 
takes to fall to 1/e times its original value, it is seen 
that the molecular relaxation times are 



T 2 = l/(ifc' + 3fc) 

r 3 =l/4£' 



(9) 



in accord with eq (6). 

Information regarding the nature of the modes of 
decay can also be obtained from the solutions. The 
constants C i2 in eq (7) show that mode 2 promotes 
the attainment of equilibrium by permitting a flow 
of dipoles from site 1 equally to sites 2, 3, and 4 (or 
vice versa) . Similarly, it is seen that mode 3 involves 
reorientations among the equivalent sites only; site 
1 is not involved. The uncertainty (which is removed 
when we place an electric field on the system) regard- 
ing the exact nature of the motions for mode 3 is a 
result of the fact that this mode is degenerate. 

2.4. Polarization and Orientational Polarizability 
Associated With the Relaxation Times 

It is desirable to know the polarization and po- 
larizability associated with each mode of decay. 



This means that we must evaluate the C# in eq (7) 
for 0^2. To do this we first calculate the new 
equilibrium number of dipoles, Ofi, in each site when 

a static electric field F is present, where F is the 
magnitude of the field. 

According to the Boltzmann distribution law we 
may write 

CfJC^CneMW/CneW&T ( 10 ) 

where C[ x is the equilibrium number of dipoles in 

site i in the presence of the field F, and C& is the 

corresponding number for site 1. The quantity fi-F 
is given by 



\fx-F) i=idF[sm £• cos 0- sin £ cos 0+ 

sin £■ sin 0- sin £ sin 0+cos £ cos 0] . 



(id 



The angles £' and 0- define the angular position of 
the ith site, while £ and define the field direction 
(fig. 2). After calculating the ratio C&/C& for r=2, 
3, and 4, and applying the condition SC/^JV where 

i 

N is the total number of dipoles in the lattice, we 
obtain 

C^Cn+lCnC^F cos SlNkT 

C^Qi-lCnCnuF cos Z/3NkT- 

C 21 »F sin Xsin£(cos0+ ^sind)/2kT 

C^ = C 2 i-^CnC 2 inF ' cos Zl'SNkT- 

C 2l »F sin Xsin £(cos 0- V3 sin d)/2kT 

C& = C 21 - 4C n C 21 nF cos HSNk T+ 
C 2 ifiF sin X sin J cos 0/fc T 





Figure 2. Arrangements of the sitesjand field direction for the tetrahedral model. 

Left diagram: Heavy arrows represent the equilibrium dipole orientations. Right diagram: Polar coordinates used 
to define the field direction. The same angles with the notation %'. and e' are used to define the equilibrium dipole 
orientations. 
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to an excellent approximation. 8 X is the tetrahedral 
angle. 

The actual number of dipoles displaced by the 
field into each site for each mode is readily obtained 
if we note, as implied in sections 2, 3, that 



A^=Cfi at t=0 
Ni=Cii at t= oo. 



(13) 



With the Cfi given in eq (12), application of these 
boundary conditions to eq (7) yields 



<7 12 =4 C u C 2l M F cos m N kT 
C 2 z= — C 21 fiFsin X sin g (cos 0+V3 sin 6)j2kT 
C zz =C 2l \x F sin X sin £ (V3 sin 0-cos 0)/2 kT 
C&=C 2l fx F sin X sin £ cos djkT 



(14) 



Now that the coefficients C# have been evaluated 
for the case of an applied electric field, the polariza- 
tion and polarizability associated with each mode 
at t=0 may be calculated. For each mode we 
obtain the components of the polarization in the 
Xj y, and z directions with 



Pz=MSMsin£cos0; 
Py=P S N t sin £ t - sin 6- 
P z =h S N t cos # 



(15) 



where the AT* (at t=0) are given by eq (7) after 
relations (14) have been inserted. To obtain P x , 
P y , and P 7 for a given mode, only the part of the 
N { involving that mode is employed. The com- 
ponents of P x , P y , and P z in the field direction are 



P(x,f)=P x sin £ cos 6 
P(v,F)=P y sin £ sin 6 
P {z ,f)=Pz cos £. 



(16) 



The average value of the polarization for all orienta- 
tions of the applied field is given by 



P$ — P{x,F) ~\~P(v,F) ~\~P(z,F) • 



(17) 



The averaging over all orientations of the applied 
field implies that the sample is poly crystalline. 

Following the above procedure we obtain for 
mode 2 



Pi- 



\ 3N 



) \3kT, 



(18a) 



8 In applying eq (10), the exponentials have been expanded and terms OQi-Fj 
k T) 2 and higher dropped. This is permissible since for ordinary field strengths 

tiF/kT«l. 



and for mode 3 



/8(7 21 \ (jtF\ 
\ 3 )\3kTf 



(18b) 



These are the polarizations referred to in eq (8). 
Using the defining relation Pp=NapF, we find the 
molecular orientational polarizabilities to be 



-W0G&K 1 "- JL 



«2 



_/8<7 2 A / m 2 \_8 
a3 ~\WjhkTr3 



iw+mu^Mf/ (18c) 



l(k> 



L1(jL\ 



+3k)J \3kTj 



(18d) 



where we have made use of the fact that Cn/N=k'/ 
(k' + 3k) and C 21 /N=k/(k' + 3k). The total polariz- 
ability, a total , is given by [8k(k' + k)/(k' + 3k) 2 } 
( f x 2 /3kT)or[Sg(l+g)/(l + 3g) 2 ](fx 2 l3kT)wheveg=k/k\ 

2.5. Characteristics of the Dielectric Relaxation 
Spectrum 

So far, we have given our results in terms of 
quantities that describe how the polarization will 
decay with time after the abrupt removal of a static 
field. It should be noted that the relaxation times 
and polarizabilities which describe the rate of decay, 
as well as the relative magnitude of each decay proc- 
ess, are expressed in terms of molecular properties. 
It is desirable to apply these results to macroscopic 
systems where the measuring field is sinusoidal so 
that we may learn something of the shape of the 
regions of loss in the dielectric relaxation spectrum. 

If the polarization of a macroscopic system decays 
after the abrupt removal of a static field with time 
according to [cf. eq (8)] 9 



P(t) = ^2P,e-^, 



(19a) 



where T$ is the macroscopic relaxation time for the 
/3th mode of decay, we may write (assuming linear 
superposition) for sinusoidal fields 



e"(co) = i:A6,coT,/(l + co 2 T^. 



(19b) 



e"(w) is the dielectric loss factor, co the angular 
frequency, and Ae# the increment of dielectric con- 
stant associated with the /3th mode of decay. Equa- 
tion (19b) reduces to the form of the Debye equa- 
tion 10 for e"(co) when there is only one relaxation 
time. For an Onsager-like internal field, it has been 
suggested by Powles n that T^r^, and we will 
employ this approximation. Also, for this field it is 
reasonable for the purposes of an approximate 
analysis to assume that Aep^K-ap where K is a con- 

9 The macroscopic relaxation time exceeds the molecular relaxation time owing 
to the tendency of the internal field to maintain itself, and the magnitude of 
the effect is proportional to the difference between the applied field E and the 
internal field F. 

10 See p. 94 of reference in footnote 3. 

» J. Q. Powles, J. Chem. Phys. 21, 633 (1953). 
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stant. 12 With these approximations we have the 
working equation/ 2,2 



€"(«)^X.;2a, W T,/(l + co 2 r£ 



(20) 



This equation, together with eq (9) and (18) can be 
used to predict the general characteristics of the 
dielectric relaxation spectrum for the tetrahedral 
model. For the present, we will restrict the discus- 
sion to the nature of a plot of e" against log co. This 
discussion will be facilitated by reference to table 1 
where the ratios a 2 /a 3 and t 2 /t s are given, together 
with values for the total orientational polarizability, 
<*totai, for various values of k/k'. It is convenient to 
discuss the results in terms of what happens as the 
ratio k/k' goes from zero to unity. 

When site 1 is relatively stable so that k/k' is 
small compared to unity, it is seen from the table 
that the relaxation times are moderately far apart, 
and that the polarizabilities are unequal but still 
of comparable magnitude. Under these conditions 
the dielectric loss peak calculated with eq (20) is 
not only broader than that predicted by the Debye 
equation with a single relaxation time, but is also 
slightly asymmetrical. Since the longer relaxation 
time, t 2 , is associated with the larger polarizability, 
a 2 , the loss peak is somewhat steeper on the low 
frequency side (fig. 3, curve A). Since the total 
polarizability is rather low when k/k' is small, that 
part of the dielectric constant resulting from orienta- 
tion of dipoles will be small. 

If we now raise k/k' to the value one-third, the 
relaxation times are somewhat closer together than 
before, and the polarizabilities are equal. The cal- 
culated loss maximum is still somewhat broader 
than one would predict from the Debye equation, 
and is also symmetrical (fig. 3, curve B). 

Table 1. Properties of the tetrahedral model 





Ratio of the 


Ratio of the 


Total polariz- 
ability atotal 
in units of 


g=k/k' 


relaxation 


polarizabilities 




limes T2/T3 


0:2/0:3 


nV3kT 





4 


2 





1/100 


3.88 


1.94 


.08 


1/10 


3.07 


1.54 


.52 


1/3 


2 


1 


.89 


1/2 


1.60 


.80 


.96 


3/4 


1.23 


.62 


.99 


1 


1 


.50 


1 



As k/k' is increased above a value of one-third, 
the relaxation times continue to converge so that 
the loss maximum becomes even narrower. The 
polarizabilities are again unequal so that the loss 
maximum again becomes asymmetrical, but owing 
to the close spacing of the relaxation times the effect 
will hardly be noticeable. It should be observed 
that this slight asymmetry is of a type opposite 
to that found for values of k/k' < 



:%. 



12 These approximations arc probably fairly accurate, especially in the case 
where the dielectric constant is low. 




Figure 3. Shape of the dielectric loss peaks for the tetrahedral 
model. 



br 



Curve A. asymmetrical and broadened peak calculated ,- <>r part); curve B, 

broadened but symmetrical peak calculated for q=Va\ curve C, narrow and sym- 
metrical peak for g = \ (single relaxation time). u> is the measuring frequency, 
and coc is the frequency where t" is a maximum. For convenience in comparing 

the shape of the curves, the maximum value for the loss has been matched at X. 



Finally, as the ratio k/k' tends to unity, the total 
polarizability reaches the full value pf/WkT, and 
only a single relaxation time T 2 = T S = %k remains. 
Thus, the Debye equation for e"(co) holds (fig. 3, 
curve C). It must be emphasized that, although 
the polarizability becomes p 2 /3kT, the value for a 
free rotator, we are dealing with hindered rotation 
inasmuch as barriers still exist between the sites. 
Hence, the observed mean relaxation time for the 
model will always considerably exceed that of a free 
rotator. 

The dependence of the shape, breadth, and magni- 
tude of the loss maximum with increase of tempera- 
ture is easilv obtained if we note that 



g= k/k'=e- v l kT . 



(21) 



It is clear from this expression that the ratio k/k' 
will increase with rising temperature 4 even if V is 
a constant. Therefore, the previous discussion of 
the changes in the dielectric relaxation spectrum 
for the case of increasing values of k/k' also holds 
for the case of rising temperature. Thus, for the 
tetrahedral model it is seen that the loss maximum 
will generally tend to become narrower with rising 
temperature. The loss peak will be steeper on the 
low frequency side if k/k' is between zero and one- 
third; i. e., if V/kT\ies between 00 and 1.1. Hence, 
for a very large range of temperature one should 
expect to find some asymmetry, and this asymmetry 
should tend to diminish as the temperature is 
increased. An analysis shows that the total polar- 
izability, and hence the dielectric constant, will rise 
up to a temperature where T^V/1.95k. Above this 
temperature the 1/Tterm in the polarizability begins 
to predominate so that the total polarizability will 
fall with further increase of temperature. 
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3. Six-position Model With One Stable Site 

The details of the model are shown in figure 4. 
The dipole in any lattice point may lie in either di- 
rection along each of the Cartesian axes, x, y, and z. 
Site 1 is unique ; all the other sites are equivalent and 
may be unstable with respect to site 1. The model 
might correspond, for instance, to the case where the 
central dipole in a body-centered cubic lattice prefers 
to point toward one face of the cube, but may also 
point toward any of the five remaining faces. 

The S-J hypothesis is assumed. A dipole can 
turn from site 1 to any of the adjacent sites (2, 3, 4, 
or 5) with a probability k where k is defined as in 
eq (1). All other reorientations to adjacent sites 
occur with a probability k' [eq (2)]. The assigned 
probabilities for the various elementary processes are 
shown in figure 4. It should be observed that the 
double jumps 1— >6, 2^4, and 3-^5 (or the reverse 
processes) are forbidden as elementary processes by 
the S-J hypothesis. 

The rate equations are 



dNJdt^-M^+k'iNz+Ns+Nt+N,) 
dN 2 ldt=kN l -^c'N 2 +k , (N 3 +N 5 +N (i ) 
dN 3l /dt=kN 1 +k'N 2 -4:k / N 3 +t(N 4 +N & ) 

dNJdt=kN 1 +k f N 3 -4,k f N,+k f (N 5 +N & ) 
dN 5 /dt=kN 1 +k f (N 2 +N i )-U / N 5 +k / N 6 
aNJdt=k'(N 2 +Nz+N 4 +N 5 )-4k'N & . 



r (22) 



(23) 



The solutions may be written in the form 

N t = On + 0^+0^+0^+ Crfi ; 

*, = «-«", 

./, _ e - (2fc + 3A/ -Q) t 

Q=^4k 2 -8kk'+5k'\ 

Substitution of the solutions CW'/* for 0^2 in the 
differential equations yields the following: 13 



Mode 2. 



^ 12=^62 — 
C22— C42 
t'32— £ ; 52- 



(24) 



This mode of decay thus involves only those dipoles 
reorienting in sites 2, 3, 4, and 5. The mode is 
active in the dielectric relaxation spectrum. 

w Mode 2 is doubly degenerate, and is handled in the same manner as modes 
3 and 4 for the tetrahedral model. 




Figure 4. Six-position model with a single stable orientation 

(site 1). 

Heavy straight arrows indicate the equilibrium dipole orientations. Light 
dashed arrows indicate elementary processes which take place with a probability 
k. All other jumps to adjacent sites (such as 2->l or 2->3) take place with a 
probability V. Jumps of the type l->6 or 5->3 are forbidden as elementary 
processes. The coordinates used to denote the field direction and orientations 
of the sites are the same as those given in figure 2 (right diagram). 



Mode 3. 



Ci3 — 063 — 

^23 ~ ^33=^43— " 



■C. 



,} 



(25) 



Again, only dipoles in the plane containing sites 2, 
3, 4, and 5 are involved. This mode of decay is in- 
active in the dielectric relaxation spectrum owing 
to a compensatory motion of the dipoles; the sum 

of the vectors C ?:3 is zero, so no polarization results. 
A similar situation arises for the single-axis rotator 
with four equivalent sites. 1 

Mode Jf. 

si si si si ^ 

( 24 — <-34 — ^44 — ^54 

Cu=-4(l-l/fl)Ci4 \ (26) 

64: =—4:C 2 4/H j 

where H=2g—l—q, q=Q/k', and g = k/k'. 
Mode 5. 

^-25— ^35 — ' 45— ^55 

Cia= -4(1-1/^)^25 Y (27) 

Ck=-4C 2: JH' 

where H' = 2g—l+q. Both modes 4 and 5 are 
complicated and involve dipoles in all of the sites. 
The polarizability associated with each mode of 
decay is readily calculated in the same manner as 
for the tetrahedral model. We summarize the 
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results below together with the relaxation times. 



a 3 =0, r 3 =l/6F; 



(28) 



a 4 



_ (/(g+3-2^[(2^-l + g)(2^+l)-2(^-l)] / M 2 \ 

AafcT/ 



(l + 5r/) 2 ^ 



T 4 =l/(2Ar+3fc , -0); 



«5= 



g(2g-3+g)[(2g-l-g)(2 g +l)-2(g-l)] / m 2 \ 



(l + 5g) 2 g 



\3kTj 



r 5 =l/(2k+3k'+Q). 



The total polarizability is given by [12g(l + 2g)/(l + 
5g) 2 ](/i 2 /3fer). Since the results are rather compli- 
cated, we resort to calculation of the polarizabilities 
and relaxation times numerically for various values 
of k/k' . The results are given in table 2. 

It is convenient to discuss the results in terms of 
what happens as k\k' goes from zero to unity. It 
is seen in table 2 that the maximum separation of 
the relaxation times occurs when k/k' tends toward 
zero. We then have t 4 /t 5 ^ 6.85. The two most 
important relaxation times, r 2 and r 4 , tend to con- 
verge as k/k' increases, with the result that the 
corresponding dielectric loss maximum will at first 
be somewhat broader than a Debye type peak, but 
becomes more narrow as k/k' increases, and finally 
assumes the simple Debye form when k—k f . The 
total polarizability always increases with increas- 
ing k/k'. 

As in the case of the tetrahedral model, the loss 
peaks will be noticeably asymmetrical for low values 
of k/k' . Owing to the fact that the longest relaxa- 
tion time is associated with the largest polarizability 
when k/k' <^0.2, the loss maximum will be steeper 
on the low frequency side in this region. At about 
k/k' = 0.2 the peak will be slightly broadened but 
nearly symmetrical since the dominant polarizabilities 
are equal. A slight asymmetry again appears when 
k/k' exceeds ~0.2, but it will not be particularly 
noticeable owing to the closeness of the relaxation 
times. It should be observed that the slight asym- 
metry that appears in this region is such that the 
loss peaks will be steeper on the high frequency 
side. The asymmetry completelv disappears when 
k = k'. 

Table 2. — Properties of the six-position model 



g=k!k' 


Relaxation times 
(in units of 1,/fc') 


Polarizabilities 
(in units of /^ArT) 


T 2 


r 4 


T. 


a 2 


«4 


a s 


"total 




1/100 
1/10 
1/5 

1/2 

1 


0. 25 

.25 
.25 
.25 
.25 

.25 


1.31 
1.25 
.88 
.66 
.39 
.25 


0. 191 
.191 
.190 
.189 
.185 
.167 



.038 
.267 
.400 
.571 
.666 



.069 
.347 
.404 
.377 
.333 



.004 
.026 
.036 
.031 





. Ill 
.640 
.840 
.979 

1 



The expression k/k'=exp(—V/kT), where Fis the 
energy difference between site 1 and each of the other 
sites, gives the dependence of k/k' on temperature. 
It follows that k/k' will increase with rising tem- 
perature, so the discussion given above in terms of 
growing values of k/k' also holds for the case of 
increasing temperature. The dielectric behavior of 
this model is qualitatively similar to that of the 
tetrahedral model. The total polarizability rises 
up to the temperature T^V/2.17k, but falls as the 
temperature is increased further. 



4. Discussion 

It is clear from the results cited above that an 
anisotropy in the crystalline field can lead to the 
existence of a set of discrete dielectric relaxation 
times for polar molecules that possess oricntational 
sites arranged in three dimensions. The separation 
of the relaxation times is a maximum when the 
crystalline field is highly anisotropic. The relaxa- 
tion times tend to merge as the temperature is 
increased. 

The following experimental predictions hold for 
both models: (1) The regions of loss in the dielectric 
relaxation spectrum will become narrower with 
increasing temperature. (2) At low temperatures 
the dielectric loss peak will be steeper on the low 
frequency side. This asymmetry will decrease as 
the temperature is raised, and practically completely 
disappear somewhat before the loss peak assumes 
the simple Debye form. At high temperatures the 
Debye equations will be a good approximation. 
(3) The dielectric constant will be small when the 
temperature is low, and will increase as the tem- 
perature is raised; at temperatures above T^V/2k, 
the dielectric constant will decrease with rising 
temperature. It should be noted that no real 
substance should be expected to show the entire 
range of behavior mentioned above since melting or 
transition phenomena may intervene. 

It should be observed that the effect of a set of 
(discrete) active relaxation times on the shape of the 
loss regions is virtually indistinguishable from the 
effect of a suitably chosen continuous "distribution" 
of relaxation times. Also, it should be noted that 
the effects that have been predicted for the case 
of rising temperature with V constant will occur in a 
more marked way if V decreases with rising tem- 
perature, as is often assumed for the case of coopera- 
tive interaction. 

The models treated here are obviously rather 
specialized, and do not necessarily correspond closely 
to any real crystal. For instance, one might expect 
the anisotropy of the crystalline field to be such 
that more than two transition probabilities are 
required to represent the system. More general 
models are being studied. 



Washington, January 5, 1955. 



363 



